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Regular arrays of electromagnetic resonators, in turn coupled coherently to individual quantum two-level 
systems, exhibit a quantum phase transition of polaritons from a superfluid phase to a Mott-insulating phase. 
The critical behavior of such a Jaynes-Cummings lattice thus resembles the physics of the Bose-Hubbard model. 
We explore this analogy by elaborating on the mean-field theory of the phase transition, and by presenting 
several useful mappings which pinpoint both similarities and differences of the two models. We show that a 
field-theory approach can be applied to prove the existence of multicritical curves analogous to the multicritical 
points of the Bose-Hubbard model, and we provide analytical expressions for the position of these curves. 
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I. INTRODUCTION 

Quantum phase transitions in interacting systems composed 
of particles of more than one species have recently produced 
a lot of interest. A host of problems in condensed matter 
physics fall into this category, prime examples being heavy 
fermions in Kondo lattices,— multi-component systems of ul- 
tracold atoms in optical lattices,**^ and the ensemble of two- 
level atoms interacting with a bosonic mode described by the 
Dicke model 

A recent addition to this list of systems is the Jaynes- 
Cummings (JC) lattice model, see Fig.[T] It describes an array 
of electromagnetic resonators, each coupled coherently to a 
single quantum mechanical two-level system (also referred to 
as "atom" or "qubit"). Recently, this JC coupling^ has been 
studied extensively in the context of quantum computation 
and quantum optics.i^ The interplay of Jaynes-Cummings in- 
teraction, observable in experiments as a photon blockade,™ 
and the transfer of photons between nearest-neighbor res- 
onators within the array render this problem nontrivial, and 
open a new route to study strongly correlated photon physics 
and quantum phase transitions of light] ^''^^-'-^-''* Similar to ul- 
tracold atom systems, the Jaynes-Cummings lattice could pro- 
vide a novel quantum simulator of condensed matter Hamil- 
tonians. 




FIG. 1: (Color online) Schematic of the Jaynes-Cummings lattice 
system, consisting of an array of electromagnetic resonators, see (a), 
with a coupling between nearest-neighbor lattice sites due to photon 
hopping. Each resonator is coherently coupled to a two-level system, 
shown in (b). 



Previous work on the JC lattice model has so far ex- 
plored several aspects of its theory, and has provided the first 
valuable hints towards possible realizations of the model in 
experiments. The early mean-field treatmentJ^ii^ has pro- 
vided evidence for a phase transition from Mott-insulating 
phases to a superfluid phase of polaritons, resembling in large 
parts the physics known from the Bose-Hubbard modelji^ 
Both numerical and analytical methods have subsequently 
been employed to extract the phase boundaries beyond the 
mean-field approximationr^ii2i22i2ii2^ and to study the criti- 
cal behavior^SiSLS^ flie spectrum of excitationsr^i^LS and the 
physics in the limit of only a few coupled cavities;^ For a 
more comprehensive review of the relatively young field, we 
refer the reader to the excellent paper by Hartmann et ali2i 

The prime motive of the present paper is to clarify the na- 
ture of the coiTespondence between Bose-Hubbard model and 
Jaynes-Cummings lattice model. After the introduction of the 
specific model in Section |II] we provide a systematic discus- 
sion of the asymptotic behavior of the JC lattice model in 
several useful limits (Section HIIl i. The knowledge of these 
asymptotics aides our discussion in the subsequent three main 
sections of our paper. The mean-field treatment is revisited in 
Section HV] where we provide analytical expressions for the 
phase boundary for general atom-photon detuning, and elab- 
orate on several issues with previous mean-field treatments. 
The correspondence beyond the mean-field level is explored 
in Sections |V] and [Vl] where we discuss several useful map- 
pings of the model and formulate the field-theoretic treatment 
of the Jaynes-Cummings lattice model. The latter approach 
allows the existence proof of multicritical curves analogous 
to the multicritical points present in the critical behavior of 
the Bose-Hubbard model. We conclude our paper with a sum- 
mary and outlook of relevant future research directions in Sec- 
tion |VlI] 



II. JAYNES-CUMMINGS LATTICE MODEL 

The Jaynes-Cummings lattice model is described by the 
Hamiltonian (h = 1) 

H = ^Hf + F''°p - nN, (1) 
i 
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where 

Hf = uja]aj + ea+aj + 9{a]aJ + a+aj) (2) 

denotes the Jaynes-Cummings Hamiltonian on site j. (A suc- 
cinct summary of this model, exactly solvable for one site and 
well-known in quantum optics, is given in Appendix lAl) The 
operator aj (a|) annihilates (creates) a photon in the resonator 
at site j, in a mode with frequency lo. Each resonator hosts a 
single two-level atom, whose transition frequency is given by 
e. Adopting the (pseudo-)spin language, the corresponding 
lowering and raising operators are denoted by the Pauli matri- 
ces crj^ . The atom-photon coupling of strength g allows for the 
coherent interconversion between photonic and atomic exci- 
tations within each resonator. We note that for most systems, 
the JC Hamiltonian (|2]i is only obtained when invoking the 
rotating wave approximation (RWA). The contributions from 
counter-rotating terms become significant when the coupling 
strength g is of the same order as the resonator or atom fre- 
quency - a fact neglected in several recent publications where 
g = e was usedJ2ia2 

The hopping of photons between nearest-neighbor cavities 
is described by 

H'""' ^~nY^{a\a,+a]a^), (3) 

(i-O) 

where the hopping strength is parameterized by k, so that the 
overall lifetime of each resonator is given by {zck)~^, be- 
ing the coordination number of the lattice. 

The chemical potential term in Eq. ([T) is in accordance with 
our treatment of the system within the grandcanonical ensem- 
ble. The chemical potential /i couples to the (conserved) num- 
ber of polaritons, iV ~ + o'/^'J )' usual plays 
the role of a Lagrange multiplier fixing the mean total number 
of polaritons on the lattice, ( A^) . By contrast to the situation 
in ultracold atom systems, it is important to emphasize that, 
for the physical realizations of the JC lattice system proposed 
so far^i the chemical potential yu is not a directly accessible 
"knob." Instead, appropriate preparation schemes have to be 
devised to access states with different mean polariton num- 
bers, see e.g. Ref.l23l 



III. ASYMPTOTIC BEHAVIOR OF THE JC LATTICE 
MODEL 

Similar to the situation encountered in the study of the 
Bose-Hubbard modelr^ the JC lattice model is not gener- 
ally amenable to an exact solution. In the following, we dis- 
cuss two particularly simple limits: the atomic limit and the 
hopping-dominated limit where Kjg <C 1 and njg ^ 1, re- 
spectively. The understanding of these limits serves as a use- 
ful building block and sets the stage for our subsequent study 
of the more intricate and interesting regimes of the JC lattice 
model. 



A. Atomic limit: n/g <^ 1 

In this limit, the photon hopping H^°'' may be treated per- 
turbatively. To leading order, the Hamiltonian decouples in 
the site index and reduces to pure Jaynes-Cummings physics, 
described by the Hamiltonian Hj'"—fj,nj and identical for each 
site j. Consequently, the ground state of the lattice system cor- 
responds to a product state | vfi )®j = | vj; )^^_^ ® \ '^)^® ■ ■ ■ 

of the ground states of each Hj'^. (To simplify notation, we 
drop the site index j for the remainder of this section.) 

The eigenstates of the JC Hamiltonian correspond to polari- 
tons - quasiparticles consisting of both photonic and atomic 
excitations. Their spectrum, shifted by the chemical potential 
according to Eq = Eq — fin and E!^^ = En± — nn, composes 
the usual "dressed-state ladder" of the JC model: £^o = and 

£;„± =nw + A/2±[(A/2)2 + ng2]i/2 (n > 1). (4) 

(For details, see Appendix lAl) Here and in the following, 
A = e — Lu denotes the detuning between the two-level atom 
and the resonator frequency. The polariton states | n± ) are 
simultaneous eigenstates of H^'^ and of the polariton num- 
ber n = cr+CT^ + fltfl. (Note that while the symbol 'W" 
refers to the polariton number on the entire lattice, "n" de- 
notes the polariton number on a single site.) Each parity dou- 
blet ± is adiabatically connected to the symmetric and anti- 
symmetric superposition of atomic and photonic excitation, 
I n± ) = (I n i ) ± I (n - 1) T ))/\/2, which form the eigen- 
states in the resonant case, A = 0. 

To specify the ground state of the system for given chemical 
potential /i, we must find the state | na ) with minimal energy, 
i.e. 

Eii^=min{E^,E^^,E^^,...}. (5) 

Immediately, one can rule out all symmetric states | n+ ) as 
ground states, since El^^ > E^_. Hence, the ground state 
is either the 0-polariton state | ) or one of the antisymmetric 
states I n— ). To further specify the ground state, first con- 
sider the Umit (w — /i) ^ g, |A|. (Both oj and g are assumed 
to be non-negative.) In this case, the first term of Eq. (HJi dom- 
inates over the second term, and the lowest energy is reached 
for the zero-excitation state | ) . By decreasing the quantity 
{u! — fi), one reaches a point where admitting an excitation 
to the system becomes energetically favorable. This happens 
precisely when Eq ~ E^_. Repeating this argument, one 
finds an entire set of such degeneracy points, determined by 
E^_ ~ E'^^^j^ . The full set of degeneracy points is given 
by 

(m - uj)/g = v/n + (A/2.g)2 -^n+l + {A/2g)\ (6) 

where n ~ 0, 1, . . ., which reduces to — lo) / g = %/n — 
VtTTT in the resonant case. Our results for the atomic limit, 
depicted in Fig. |2ja), are consistent with results previously 
presented by Greentree and coauthors.— 

As shown in Fig. [2a), the spacing between degeneracy 
points becomes closer as (y^t — lo) is increased towards zero. 
More precisely, by inspection of Eq. ^ one finds that for 








FIG. 2: (Color online) Properties of the JC lattice model in the resonant case, A = 0. (a) Ground states of the JC lattice system in the atomic 
limit, K = 0, as a function of detuning A. The ground state of the system is given by a product state of Jaynes-Cummings eigenstates on each 
lattice site j . Depending on the chemical potential, the system assumes either the state | ) ®^ or one of the antisymmetric states | n— ) ®^ . 
Degeneracies between | n— ) and | (n + 1)— ) mark the onset of superfluidity, occurring for finite photon hopping, k > 0. The onset points, 
for zero detuning located at (/i — ^)/g = \/n — \/n + 1, become dense as [i approaches lo. For ^ > ui, the system becomes unstable. 
Degeneracies occur at the same chemical potential for negative and positive detunings A (solid curves), except for the lowest degeneracy 
between | and | 1— )®3 where the A > case is given by the dashed curve, (b) Mean-field phase diagram of the resonant JC lattice 
system as a function of the effective chemical potential, ^ ~ uj, and photon hopping strength k. The color scale shows the magnitude of the 
order parameter ^ = (a). The value of i/; reveals the Mott-insulating phases (denoted "Ml") with ^ = and fixed number n = 0, 1, ... of 
polaritons per site, and the superfluid phase ("SF") with > Q. The phase boundary can be obtained analytically, cf. Eq. l l20t . and is marked 
by a black curve. For sufficiently large photon hopping strength, the system becomes unstable with respect to addition of polaritons. A crude 
estimate of the onset of instability is given by — uj) = ZcK, depicted by the white dashed curve. In the hatched region close to /i — tj = 0, 
numerical results are unreliable when using a fixed cutoff for the maximum photon number, (c) Improved numerical results near fi — uj = 
[range of /j, — u and k as marked by the rectangle in panel (b)] can be obtained by employing a "sliding" truncation (see text) centered at the 
photon occupation number obtained in the atomic limit. 



n oo, (ii—uj) tends to zero, independent of the detuning A. 
In other words, the point (/i — lu) = Q where the chemical po- 
tential is tuned to exactly compensate the resonator frequency, 
plays the role of an accumulation point. In SectionlTVlbelow. 
we will demonstrate that the existence of this accumulation 
point necessitates extra care in the numerical analysis when 
approaching — oj) = 0. 



This leaves the question of what happens for (/i — cj) > 0. 
In this case, the system becomes unstable: adding more po- 
laritons to the system always lowers the total energy. This be- 
havior, if not accompanied by an additional mechanism that 
would ultimately limit the polariton number, is unphysical. 
Hence, we restrict ourselves to values of /i e [— oo, uj] in the 
following. We will encounter a similar instability in the fol- 
lowing section describing the limit of large photon hopping. 



B. Hopping-dominatedlimit: k/p ^ 1 

For large photon hopping, n/g ^ 1, the Hamiltonian iJ^°P 
overwhelms the atom-photon coupling, and the latter may be 
treated perturbatively. As a crude approximation, we con- 
sider the order g'-\ i.e. atom-photon coupling is dropped and 
the atomic and photonic systems decouple completely. The 
ground state energy then consists of the atomic contribution 
(all atoms occupying their respective ground states) and the 
photonic contribution from the bosonic tight-binding model 

» ihj) 

As usual, the tight-binding model may be diagonalized in 
terms of single-particle Bloch waves. For the relevant cases 
of a 2d cubic lattice and a 2d honeycomb lattice, the resulting 
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energy dispersions are given by 

e(k) = {lo — fj.) — 2k ^ cos(fcia) 



for the cubic lattice,— and 



e-i-(k) — {uj — fi) ± K 



(8) 



(9) 



for the honeycomb lattice (a denotes the lattice constant).— In 
the energy dispersion of the honeycomb lattice, the two signs 
refer to the lower and upper (tt and tt*) band. Independent of 
the specific lattice type, the bosonic ground state is obtained 
by iV-fold occupation of the k = (tt) state with correspond- 
ing energy 



N{lo - m) - NZcK. 



(10) 



This result reveals the presence of an instability at large pho- 
ton hopping strength: for > — /i the ground state en- 
ergy becomes negative and can be made arbitrarily small by 
increasing the polariton number TV. (We note that a similar 
instability formally exists for the Bose-Hubbard model in the 
region of negative chemical potential.) This instability is al- 
ready present in a mean-field treatment of the JC lattice sys- 
tem, where its presence has apparently been overlooked in 
Refs.[iland|23 



C. Effective Hubbard- (7 in the dispersive regime? 

A prime motive of the present paper is the analysis of differ- 
ences and similarities between the Bose-Hubbard model and 
the Jaynes-Cummings lattice model. As we shall demonstrate, 
based on mean-field theory, exact and approximate mappings, 
as well as field theoretical methods, there exist strong par- 
allels between the two models. We find, however, that the 
direct comparison of the two models in terms of an "effec- 
tive Hubbard-C/" is generally not appropriate. Correcting for 
a missing factor of 1 /n in Ref. [isl one may attempt to define 
such an effective Hubbard-JJ via the relation 



C/n± = (i^f'^+Di-i^.'^i+Ml/^- 



As argued by Greentree et al., this effective U approaches 
zero in the limit of large polariton numbers n per resonator 
- consistent with the predominance of the superfluid phase in 
this regime. By contrast, we find that the system does not 
approach a constant Hubbard-J7 in the large detuning limit 
(as stated in Ref. [Tsl) when using the correct expression ( fTTI ). 
To see this, we consider the dispersive regime g/A <C 1 
with polariton numbers below the critical photon number, 
n <C nc = A^/4g^. In this case, the dressed-state ladder 
decomposes into two nearly harmonic ladders shifted by the 
detuning A, 



1^ 
2 



(12) 



i.e. photons and atoms are effectively decoupled. This sit- 
uation is not correctly captured by a Bose-Hubbard model, 
and indeed we find that the effective interaction J7„± — cu/n 
would need to be strongly polariton-number dependent in this 
case. 



IV. MEAN-FIELD THEORY 

The mean-field theory (MFT) of the Jaynes-Cummings lat- 
tice model has previously been discussed by Greentree et al.— 
It reveals a second-order phase transition from a superfluid 
phase of polaritons to Mott-insulating phases. We revisit this 
treatment, pointing out several important issues not addressed 
in RefE 

Starting from the JC lattice Hamiltonian, Eqs. ([T])-©, the 
MFT is constructed by a decoupling of the lattice sites. We 
apply the general prescription of replacing an interaction term 
involving two operators A and B by its mean-field expression, 
AB {A)B + A{B) - (A) (B), to the photon-hopping term 
in the Jaynes-Cummings lattice model. This yields 

(ij) I jGnn(j) 

i jenn(i) 

where nn(z) denotes the set of sites that are nearest neighbors 
to site i. Translation invariance implies that the expectation 
values are actually site-independent. Introducing the order pa- 
rameter as i/' = ZcK(ai),— the mean-field Hamiltonian reads 
iJ™*' = J2j hf with 

hf = i(e - pt)(T| + (w - l-L)a]aj + g{a]aj + a+ a^) 

-{a,r +a]^) + — \i^? ■ (14) 

Here, Zc is the coordination number of the lattice. It is useful 
to note that the dimension and precise geometry of the spatial 
lattice underlying the JC lattice model enter the mean-field 
theory only through Zc- Moreover, a change of the coordina- 
tion number is equivalent to a rescaling of the photon hopping 
strength. Specifically, within MFT the substitution Zc rzc 
has the same effect as the rescaling n ^ rn while keeping Zc 
fixed. 

Throughout this paper, we focus on the quantum phase tran- 
sition of the JC lattice model at zero temperature. In this 
case, the order parameter is determined by minimization of 
the ground state energy -E'o(V') of the mean-field Hamilto- 
nian h™^. Figure |2|b) presents numerical results from this 
analysis. Consistent with numerical results previously pub- 
lished by Greentree et al.r^ our phase diagram shows lobes 
where the order parameter t/j vanishes. These lobes, which 
are the analogs of the Mott-insulator (MI) phases in the Bose- 
Hubbard model,— border on a superfluid phase with nonzero 
order parameter, \'tl>\ > 0. The onset of superfluidity is con- 
tinuous, in line with the assumption of a second-order phase 
transition. 
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In addition to the Mott-insulating and superfluid phases dis- 
cussed by Greentree et al., our results also reveal the onset 
of the polariton-number instability expected for the hopping- 
dominated regime, see Section UlI Bl above. The unstable re- 
gion can be identified by the darkest coloring in Fig.|2b). In 
the numerical analysis of the mean-field equations, the insta- 
bility manifests itself in solutions with excessively large val- 
ues of the order parameter More specifically, in the un- 
stable region, the value of the order parameter is entirely de- 
termined by the photon cutoff, i.e. the truncation of the pho- 
ton Hilbert space necessary in the numerical treatment. When 
truncating this Hilbert space at some maximum photon num- 
ber n,nax, one can confirm that the order parameter obeys the 

1/2 

relation ij; ~ rimax in the unstable region. This instability ap- 
pears to be present in the numerical results obtained by Green- 
tree et al.,— and very likely also in qualitatively similar re- 
sults obtained by Lei and Lee for the Dicke-Bose-Hubbard 
model;^ In both references, the instability is not mentioned, 
and it has not been discussed elsewhere to the best of our 
knowledge. A rough estimate for the border of the instabil- 
ity region may be obtained from the hopping-dominated limit 
by setting the ground state energy, Eq. (fTOl i. to zero. This 
leads to the relation (/i — cj) = —ZcK,, shown as the dashed 
curve in Fig.|2b), which is in reasonable agreement with the 
numerical results. 

There exists a second region in the phase diagram which is 
highly sensitive to the specific photon truncation scheme. This 
region close to /i — a; ~ 0, marked by hatching in Fig. I2b), 
appears to suggest the emergence of a large Mott-insulating 
phase reaching from /i — w ~ —0.1 all the way into the 
ji — Lo > Q region. This featu re, p resent also in all mean-field 
phase diagrams of References! 1 5l andlz?! turns out to be a pure 
artefact of the truncation of Hilbert space at some maximum 
photon number. Indeed, its presence would be rather diffi- 
cult to reconcile with the atomic-limit prediction of an infinite 
number of MI phases with /i — u = being an accumulation 
point and with the region — > being unstable. 

To illustrate this point, we have employed an alternative 
"sliding" truncation scheme based on the number of excita- 
tions in the atomic limit. Specifically, if the ground state in 
the atomic limit is given by | n~ )®\ we take into account a 
total of 2 A„ photon Fock states, with numbers n' within 



max{0, n — A„} < n' < 



(15) 



The results from this sliding truncation scheme, presented in 
Fig.[2jc), uncover that results from the maximum photon num- 
ber truncation scheme are unreliable in the hatched region. As 
expected from the atomic limit, the sliding truncation scheme 
shows no large Mott-insulating region close to fi — ui ~ 0, but 
only the sequence of ever smaller MI phases approaching this 
chemical potential value. 

The phase boundary between MI and superfluid phases can 
be determined in a way similar to the procedure for the Bose- 
Hubbard model. '^■^'^ In the critical region, the superfluid order 
parameter ip is small and the terms in /i"^ involving ip may 
be treated perturbatively. This leads to an expansion of the 



ground state energy in powers of ■tp. 
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EoW = + r IV-r + 2^ IV^r + Om°). (16) 

Note that the invariance of h^^ under the global gauge trans- 
formation 



-i(f) 



— —la. 



(17) 



implies that all odd-order terms in the expansion ( fT6b vanish. 

Eq. (fT6l ) represents the standard situation of a quadratic 
plus quartic potential, ubiquitous in the study of mean-field 
phase transitions. The resulting mean-field critical exponents 
are /? = 1/2 for the critical exponent of the order parameter, 
ip ^ \jP , and a = for the exponent of the compressibility, 
defined by Kg ^ j22, Here, /2 = [ji — jjLc) / jic measures the 
relative chemical potential difference from the critical value. 
Other critical exponents can be established via scaling and hy- 
perscaling relations, and we discuss the role of the dynamical 
critical exponent z in Section lVTl 

As usual, the phase boundary is specified by the condition 
r = (assuming that u > 0). The quadratic expansion coeffi- 
cient r can be expressed as r = i?„ + (zck)^^ where i?„ can 
be obtained from second-order perturbation theory in the pho- 
ton hopping. In the resonant case (A = 0), the contributions 
are given by 



Rn=o — y2 — ; — , 



(18) 



(V" + oiVi^ — 1)^ 



— {^i — uj) — {y/n — OL\Jn — l)g 



uj) ~ (^/n — a\/n + l)g 



(19) 



where n = 0,1,... enumerates the sequence of Mott- 
insulating lobes. For nonzero detunings, the lengthier expres- 
sions for the quantities i?„ are relegated to Appendix |B] The 
full phase boundary is then obtained analytically as the set of 
curves defined by 



(20) 



z B ' 



shown for the resonant case as black curves in Fig.|2jb), and 
for finite detunings in Fig.|3] We will return to the question of 
phase boundaries in Section |Vll where the presence of multi- 
critical hnes will be discussed. 



V. USEFUL MAPPINGS FOR THE JC LATTICE MODEL 



The qualitative similarity in the phase diagrams of the 
Jaynes-Cummings lattice model and the Bose-Hubbard model 
has been one of the main driving forces in recent studies of 
the JC lattice problem. Both models display a quantum phase 
transition from a superfluid phase to Mott-insulating phases. 
Yet, as discussed in Section ITlI CI it is in general not possible 
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proper fermionic description. Here, we instead employ the 
simple equivalence between spin-1/2 systems and hardcore 
bosons, valid for lattices in arbitrary dimensioni^ Accord- 
ing to this equivalence, the energy spectrum of the spin- 
Hamiltonian H{cf'^ ,aj) is identical to the set of finite 
eigenenergies of the hardcore boson model 



U 



(23) 



where the first term is obtained from the Hamiltonian H by 
replacing all spin operators with boson operators according to 
(j^ ^ bjj, (tJ ^ bj. In our case, the resulting two-component 
Bose-Hubbard model reads 



H' = Y^ {uj- fi)a]aj + (e - fi)b]bj + -^mj {ubj - 1) 



g{a%+h\aj) 



(24) 



FIG. 3: (Color online) Mean-field phase boundary as a function of 
effective chemical potential (/i — uj) and atom-resonator detuning 
A. The cross-section shows the Mott-insulating lobes at zero de- 
tuning (in orange/light gray). With exception of the | 0) lobe, all 
lobes become narrower when increasing the detuning A. Arrows 
and curves in red/dark gray mark the positions of multicritical lines 
where the system switches from the universality class of the generic 
superfluid-Mott insulator transition to the universality class of the 
soft-spin (d + 1) -dimensional XY model. 

to view the JC lattice model as a simple Bose-Hubbard model 
with an effective Hubbard interaction U. Instead, the compos- 
ite system of interacting photons and atoms can be mapped to 
a two-component Bose-Hubbard model or a general polariton 
model which reduces to the spin-1/2 quantum-XX model in 
the vicinity of the degeneracy points between adjacent Mott 
lobes. We present the corresponding mappings and discuss 
their implications in the following three subsections. 



where Ubj = bjbj is the number operator for the b bosons on 
site j, and Uoo is a shorthand for the large-[/ limit, cf. Eq. 
( l23b . Such two-component Bose-Hubbard models are also of 
interest in systems of ultracold atoms in optical trapsi^*^ In 
that context, the two components correspond either to differ- 
ent species of atoms, or to two internal states of one atom 
species. We note that for cold atoms, the two-component 
Bose-Hubbard model has mainly been studied in the limit 
of negligible inter-component conversion, i.e. in our notation 
g = 0. The relevance of atom-photon coupling in the JC 
lattice model serves as a motivation to explore the physics 
of the two-component Bose-Hubbard model in the regime of 
g ^ and large repulsion Uoo for one of the two compo- 
nents. We note that the hopping-dominated regime is eas- 
ily obtained from ( l24l i by treating inter-component conver- 
sion perturbatively, such that atoms and photons effectively 
decouple. This results in a superfluid of photons only weakly 
dressed by atomic excitations. 



A. Mapping to a two-component Bose-Hubbard model 

As the first useful reformulation of the JC lattice model 
we explore its mapping to a two-component Bose-Hubbard 
model. For this mapping, the annihilation and creation oper- 
ators for spin waves are identified as Cj ~ a~ , = a^' . As 
usual, these operators are neither bosonic nor fermionic: while 
on each individual site the cj obey anticommutation relations, 
CjS from different sites commute, 

{c„c,} = 0, {c,,c]} = l (21) 
[c, Cj] = [c,, c]] = for i 7^ j. (22) 

This deficiency may be cured by the Jordan- Wigner 
transformation^ and its generalizations^^ii^ which yield a 



B. Polariton mapping 

The polariton mapping of the JC lattice model was pre- 
viously discussed by Angelakis et ali^ In that work, the 
mapping was considered for the situation close to resonance, 
A ~ 0. Here, we extend the mapping to the general case of 
nonzero detuning. 

The idea of the polariton mapping is to diagonalize the lo- 
cal Jaynes-Cummings Hamiltonian Hj'^ for each site j, and 
to rewrite the remaining photon hopping H^°^ in the form of 
polariton hopping between nearest neighbor sites. Hence, the 
starting point for this mapping is the atomic limit, where the 
system's eigenstates are dressed-state polaritons of the form 

|n,±)=a„±|ni)+5„±|(n-l)t). (25) 
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For n > 1 the coefficients read 

_ J sin^n, a = + 

COS U„, a = — 




a 
a 



(26) 



\n+) 








[ \n-) 




[ an- 





i("-i)T>;' 



wliile for tlie nondegenerate ground state n = we introduce 
the conventions | 0— ) = | j ) and | 0+ ) = 0, i.e. ao_ = 1, 
ao+ = bo± = 0. (For the definition of the mixing angles 6'„ 
see Appendix IaI) With this, the Hnear relations between the 
basis states may be written compactly as 



(27) 



(28) 

As required, the JC part of the Hamiltonian is now diagonal, 

oo 

'Hf~l^N = Y^Y.ll K^cPlncP^.n. ■ (29) 



with inverse 



ni) ^ 






On- A 


{\n+)\ 


(n-i)T> J 


1 = 


V 


bn- J 


[\n-) 



=0 a=± 



Here, we have introduced the polariton operators in the same 
way as defined in Ref . |23l namely 



(30) 



and their Hermitian conjugates. The operator Pj maps the 
polariton ground state | 0— )j on site j to another polariton 
state I na )j on the same site. Note that these operators do not 
satisfy the canonical commutation relations of creation and 
annihilation operators. 

To rewrite the hopping contribution in the polariton basis, 
we observe that the effect of a photon annihilation operator 
can be expressed in the polariton basis via 



\n- 
\n- 



I ng ) 

\n-l,e) 



^n,-\ — h ^n,-\ — 
^n, — + ^n, 



\n-l, + ) 



(31) 
(32) 



where the conversion amplitudes are 



tn,±+ = \/nan±an-i^+ + \Jn- 16„±a„_i,_, (33) 
tn,±- = Vnan±bn-i,+ + \/n- 16„±6„_i _. (34) 

As a result, the representation of the photon annihilation op- 
erator in the polariton basis reads 



a-j — E E in,aQ'^'2(„_i)Q/-P;,nQ, (35) 
n— 1 a,CK' — lb 

meaning that aj diminishes the polariton number n on site j 
by one, and may or may not alter the polariton type a = ±. 



In full, the JC lattice Hamiltonian can thus be rewritten as 

(36) 



j n— a — di 



-'^E E E E *n,aa'tn',/3fl' 
j j'gnn(j) ri,n' = l Q,a',/3,/3' 

^ P],na Pj,(n-l)a' Pf , („' _ i)/3' Pj' ,n'p, 

where the polariton hopping indeed transfers polaritons be- 
tween nearest-neighbor sites, and additionally may change the 
polariton type (±), of one or both sites involved. 

It is instructive to compare this form of the JC lattice Hamil- 
tonian to the Bose-Hubbard model. This comparison is facili- 
tated by rewriting the Bose-Hubbard model H = Hbh + Hhop 
in the corresponding form, 

^ oo 
HbH = E + ot^"^ ("^ - 1) = E E EnPlnP,,n, 



j «=0 



(37) 
(38) 



i?hop = -KE E 

j j'eim(j) 

= -«E E Y.^^4nPl^n-lPyn'-lPr,n', 
j i'Gnn(j) ri,Ti' 

with En = ne + Un{n — 1) /2. Inspection thus shows that the 
structures of the JC lattice Hamiltonian and the Bose-Hubbard 
Hamiltonian share several crucial properties. They both in- 
volve an onsite contribution and a hopping term which trans- 
fers excitations between nearest neighbor sites. The two sys- 
tem differ, however, in (i) the specific onsite eigenenergies, 

(ii) the scaling of the hopping matrix elements with n, and 

(iii) in the case of the JC lattice model, the existence of two 
excitation species and the possibility of interconversion. 



C. Mapping to the spin-1/2 quantum-XX model 

Angelakis et al. have pointed out that the interconversion 
between the two polariton types may be neglected within 
RWA.-^^ Specifically, their argument is based on the restriction 
to the subspace consisting of ground state and one-excitation 
manifold, and results in a mapping to the XX modelj^ We 
review this statement, specify its range of validity and gener- 
alize it to include higher-excitation manifolds. 

Consider the regime close to a degeneracy point between 
two Mott-insulating lobes (polariton occupation numbers n — 
1 and n) in the atomic limit, i.e. n/g <C 1. The effective 
model for low-energy states in this regime treats the two rel- 
evant product states | {n — 1) — and | n— )'*^ . Given the 
conditions 



K < min{|i?'' 



(n+l)- 



IE'' 



- E^ 



and 



(n-l) + 



E'; 

(n- 



W (39) 



(40) 
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the hopping of photons will not lead to significant occupation 
of the symmetric Jaynes-Cummings states (higher in energy) 
or to occupation of adjacent Mott-insulating states. Thus, we 
may truncate the Hilbert space and drop inter-species conver- 
sion, which results in 

^y~! ^{n+m)~-^l{n+m)~^3,{n+m)- (41) 

j m— — 1,0 

j i'enn(j) 

Note that the onsite energies E^_^_^^_, and the hopping 
amplitude ntn = Kt„ .-- depend on the detuning A between 
resonator frequency and qubit frequency. Changing the detun- 
ing can thus be used to effectively alter the hopping strength, 
which provides a convenient experimental "knob" to tune the 
system across the phase transitioni^^ 

Within the truncated Hilbert space, the effective Hamilto- 
nian iJ^*' can be rewritten as a spin-lattice XX model. To see 
this, one identifies the polariton operators with spin- 1/2 oper- 
ators in the following way: 



-1)- 



(42) 



With this definition, aj and its Hermitian conjugate obey 
the standard commutation and anticommutation rules for the 
Pauli lowering and raising matrices. Specifically, we have 



l"->,(("-i)-L-,l("-i)->,-(^-Lv 



= <5,,,1, (43) 
where we have used the completeness relation 

1 = I >, ( n- I,. + \{n- 1)- ( {n - 1)- \^ (44) 

valid in the truncated Hilbert space of site j. Similarly, we 
confirm that 



{a+,a+} = 0, W+^aJ} 
Together with the usual definitions 



a. 



(45) 



1, 

(46) 

Equations ( |43l l and ( |45] l are sufficient to reproduce the com- 
mon algebra of Pauli spin matrices, i.e. for a,b,c G {x, y, z} 



(47) 



The resulting XX model describing the physics close to the 
degeneracy points of the JC lattice model is given by 



(n-l)- 



(48) 



- -Ktl V 

{],]') 



This effective model underlines the similarity between the JC 
lattice model and the Bose-Hubbard model in the hardcore 
boson limit U / n oo, for which a similar XX model is ob- 
tained at the degeneracy points of adjacent Mott lobes3 



VI. EFFECTIVE FIELD THEORY AND EXISTENCE OF 
MULTICRITICAL POINTS 

We now turn to the field-theoretic treatment of the JC lattice 
model. Employing imaginary time coherent-state functional 
integration, the partition function of the JC lattice system may 
be expressed as 



Z = J llVa*iT)Va,iT)VN,iT)SiN^^ ~ 1) 



where the action is given by 



(49) 



(50) 



S[a*,a„N,] =Sb + £ | E 
j 

Here, H^'^ and iJ^°P are now the Hamilton functions obtained 
from the Hamiltonians, Eqs. (|2]i and (|3]l, by replacing all op- 
erators with their associated fields aj, a*, and Bloch vectors 

N, = (iV,-,,iV,,„iV,,,),i.e. 

flj aj{T), a*{T), <jf Nj^a, (51) 

where a = x,y, z. For convenience, we have absorbed the 
chemical potential term into the JC Hamiltonian, Hj^ 

— i^LUj. Further, the Berry phase contribution to the action 
from the spins is given by 

^^=E/"'^^(N,(r)|^|N,(r)). (52) 

(For a comprehensive exposition of coherent state functional 
integration of bosons and spins, see e.g. Refs.l28l and[34l) 

Analogous to the treatment of the Bose-Hubbard 
model^ia^ we next use a Hubbard-Stratonovich trans- 
formation to decouple the hopping term. 



exp 



/3 



3,j' 



(53) 



j n -Dr, {r)^i (r) exp - f dr ^>7,'^.- 

3 L jj, 



X exp 



Here, the nearest-neighbor hopping has been encoded in the 
matrix Kjj', yielding the photon hopping rate k if j and j' are 
nearest-neighbor sites, and zero otherwise. The introduction 
of the auxiliary fields ipj thus results in an effective action that 
is completely local. In full, the partition function now reads 

Z=j X{Vi^*{T)Vi,,{T)Va*{T)Va,{T)Vn,{T)5m - 1) 



X exp I — S" ['(/'!, V'ij Q^j I 



(54) 
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with an action 5" that includes the auxiliary fields 

S' = Sb+ E«.V + E^K.'^/ (55) 



In the next step, we integrate out the fields a*, aj, and Nj. 
For the discussion of the quantum phase transition, we are 
primarily interested in an effective field theory for the auxil- 
iary fields ip* and ipj, which should hold for "small fields" ip. 
(The fields -ipj are proportional to the order parameter {aj), 
and hence are small in the critical region.) Thus, we may ex- 
pand in the fields ijjj and their temporal and spatial gradients 
to obtain the general effective action 

S,ii[r.i^]^j^dT J d''xC,ffirix,T),^Pix,T)), (56) 
with Lagrangian density 



■r\ip\" 



dr 



dip 



dr 



•XalWr (57) 



Next, we wish to obtain explicit expressions for the coeffi- 
cients in this expansion. As in the Bose-Hubbard modelr^ia^ 
the coefficients r and u can be related to the mean-field coef- 
ficients r and u, Eq. ( fT6l ). To see this, consider the situation of 
a spatially and temporally constant field ipjir) = tjj. In this 
case, we identify 



-Serr(^) 



(58) 



where the effective action ( |56] l reduces to 

S^ffW = 0V (i^o + r + ^ IV.!"^ + •••). (59) 

Here, V denotes the volume of the system. The substitution of 
the global field ^ into the action S', Eq. (l55T l, is facilitated by 
the relation 'J2ij i^i/ = A/"/ (kz^,) where M is the total number 
of sites. By inspection, one can verify that Z{'tj}) is the parti- 
tion function corresponding to the mean-field Hamiltonian 

3 

+ g{a]<jj + a+a,) - {a,r + a>) + — l^l' • (60) 

In the limit of low temperature, /? ^ oo, and close to the 
phase boundary (i.e. small ip), we thus identify 



where Eo{ip) is the ground-state energy of the mean-field 
Hamiltonian H^f. As a result, we conclude 



u, 



(62) 



where v = Vj M is the volume per lattice site. 

To obtain the coefficients K\ and A'2, we generalize the 
time-dependent U(l) symmetry encountered for the Bose- 
Hubbard model-- to the JC lattice model. In the case of the 
latter model, the symmetry involves both the photon fields as 
well as the spins: 



(63) 











/ cos (/3 


sinyj 




-> — sin if) 


cos ip 




\ 





i'3 - 








(cj - /i) - 


dtp 



This symmetry transformation leaves the action S', Eq. ( |55] l 
invariant^ The only two terms for which invariance may not 
be immediately obvious are the atom-photon coupling term 
and the spin Berry phase term Sb- The verification of invari- 
ance for these two terms is given in Appendixicl 

With the invariance intact, we now explore its consequences 
for the expansion coefficients in Eq. ( |57] |. Specifically, we 
consider the transformation in the limit of a small phase 
(P{t) <C 1, and plug the transformed fields tp*, ip, and the 
transformed lu — ji into the effective action S'eff[V'*, Tak- 
ing into account that the coefficients f , u, and Ki may depend 
on a; — yu, we Taylor-expand in powers of the phase ip and its 
time derivatives. Requiring that S'eff be invariant, the leading 
order in this expansion yields the important relations 



dr 



d{uj — ^) ' 



2 d{LO - 



(64) 



(61) 



see AppendixiDjfor a more detailed derivation. 

Whenever the coefficient Ki vanishes, the phase transi- 
tion of the JC lattice system changes its universality class. 
While the case of Ki ^ is associated with the generic 
superfluid-Mott insulator transition with dynamical critical 
exponent z = 2, points where Ki = Q lead to a phase tran- 
sition within the universality class of the soft-spin {d + 1)- 
dimensional XY model. The physics of these multicritical 
points is in complete analogy to the corresponding physics 
of the Bose-Hubbard modelfii^S 

One difference worth noting is the presence of an additional 
parameter in the case of the JC lattice system. This parame- 
ter, the energy scale for the atoms e, is to be held constant 
when calculating the partial derivative in Eq. ( |64l i. Due to the 
presence of this additional parameter, the phase boundary is a 
two-dimensional surface in the space spanned by the parame- 
ters Lo, £, and K, and the condition if 1 = defines curves on 
the phase boundary. Consequently, the field theory predicts 
a multicritical curve, as opposed to the isolated multicritical 
points of the Bose-Hubbard model. 
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FIG. 4: (Color online) Change of variables shifts multicritical points 
away from lobe tips. While multicritical points appear at lobe tips 
when plotting the phase boundary as a function of {u — /i) and n for 
constant atomic energy e, their position is shifted when plotting the 
phase boundary for constant detuning A, see Eq. ( I65t . 



These curves are shown in Fig. [3j where the axes have 
been chosen as resonator frequency uj and detuning A. As 
mentioned by Zhao et al.,^" the transformation of variables 
{a;,£} — > {tj, A} shifts the position of multicritical points 
away from the lobe tips, when considering cuts for constant 
detuning A, compare Fig.|4] Specifically, the change of vari- 
ables implies the transformation of the partial uj derivative ac- 
cording to 



(65) 



where subscripts to the partial derivatives denote which vari- 
able to hold constant. 

We thus predict multicritical curves, whose position is com- 
pletely specified through the conditions given by Eq. ^'^i^ 



d 




d 




d 




€ 




A 


OA 



d_ 

duj 



Rn = 0. 



(66) 



For each Mott lobe with n = 1, 2, . . ., there is one such mul- 
ticritical curve, along which the universality class changes in 
a way similar to the Bose-Hubbard model. It is interesting to 
note that the numerical evidence for these multicritical curves 
is currently discussed controversially. Specifically, Zhao et al. 
argue for the complete absence of multicritical points based 
on their quantum Monte Carlo simulations.-' By contrast, 
Schmidt and Blatter have recently presented new evidence for 
the presence of such multicritical points/s^ which is consistent 
with our findings from field theory. While the exact position 
and shape of the phase boundary may acquire corrections be- 
yond the mean-field results presented in Eqs. (fT8]l-(|20]i. we 
emphasize that the field-theoretical analysis is general and ap- 
plies to arbitrary spatial dimensions of the system. Specif- 
ically, as long as the phase boundary remains differentiable 
(i.e. no kinks) it predicts the existence of multicritical lines 
independent of the spatial dimension of the lattice. 




sc ( ubit 



sc resonator 



FIG. 5: (Color online) Possible realization of the JC lattice model, 
building upon circuit quantum electrodynamics: superconducting 
microwave resonators can be patterned onto a chip and coupled to 
superconducting qubits (shown in the zoom-in). 



VII. CONCLUSIONS AND OUTLOOK 

In conclusion, we have presented a thorough analysis of 
the quantum phase transition predicted for polaritons in the 
Jaynes-Cummings lattice. We have revisited the mean-field 
phase diagram, clarified the existence of an unstable re- 
gion, and derived analytical expressions for the full two- 
dimensional surface that constitutes the phase boundary. The 
qualitative similarities between the Jaynes-Cummings lattice 
and the Bose-Hubbard model have been elaborated and dis- 
cussed in the context of exact mappings to a two-component 
Bose-Hubbard model and a generalized polariton model. Fi- 
nally, we have presented a field theoretic approach for the 
Jaynes-Cummings lattice, whose gradient expansion in the 
critical region underlines the analogy with the Bose-Hubbard 
model. We find that the Jaynes-Cummings lattice falls in the 
same universality class, and have proven the existence of mul- 
ticritical curves which parallel the presence of multicritical 
points in the Bose-Hubbard model. 

A number of potential candidates have been suggested for 
an actual realization of the Jaynes-Cummings lattice physics, 
ranging from arrays of photonic band-gap cavities to realiza- 
tions in circuit QED; see the review by Hartmann et al. and 
references therein.-"* The realization in circuit QED, see Fig. 
|5] is particularly interesting. In this case, the fabrication of 
the basic building block is well established and shows the de- 
sired Jaynes-Cummings physics.i^^ The step to fabricating 
medium-size arrays {M ^ 10^-10^) should not pose fun- 
damental difficulties, and first ideas how to initialize such a 
system have been suggested.-- Namely, the system could be 
prepared within a Mott-insulating phase by using a global ex- 
ternal microwave signal. An off-resonant ac Stark tone could 
then be used to tune the system across the phase boundary by 
modifying the detuning A, and would also allow for the study 
of the system at multicritical points. By contrast, both k as 
well as Lo would typically be fixed by the fabrication param- 
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eters. The order parameter ip = (a) could be accessed by a 
homodyne measurement of the voltage field in the resonators. 

In the future, these ideas will have to be extended to specify 
actual experimental protocols. As pointed out previously^^ 
an analysis of the Jaynes-Cummings lattice when subject to 
unavoidable decoherence mechanisms will be crucial to aid 
the experimental attempt. An additional layer of interest- 
ing physics will arise from disorder Specifically, in the cir- 
cuit QED architecture such disorder is likely to occur for the 
atoms' frequencies due to the limited precision in Josephson 
junction fabrication. We are confident that a concerted ef- 
fort of theory and experiment can lead to the realization and 
further study of Jaynes-Cummings lattice physics in the near 
future, and note that such a realization would share many of 
the interesting aspects and benefits of the thriving field of ul- 
tracold atom physics. 



eigenstates corresponding to these eigenenergies are | | ) for 
the ground state and 

I n- ) = sin6i„ I n J, ) + cos6'„ I (n - 1) t ) , (A4) 
I n+ ) = cos0„ I ?i i ) - sin0„ I (n - 1) t ) . (A5) 

The mixing angle 9n depends on detuning, coupling strength, 
and polariton number: 

0„ = iarctan(^^). (A6) 



APPENDIX B: PERTURB ATI VE SHIFTS FOR A / 
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In the critical region and within mean-field theory, the en- 
ergy shifts Rn IV'I^ due to photon hopping can be calculated 
via perturbation theory. The resulting expressions for vanish- 
ing detuning were given in Section HV] Here, we provide the 
slightly lengthier but general expressions: 



APPENDIX A: SOLUTION OF THE JC MODEL 

For completeness, we compile the most important results 
for the Jaynes-Cummings model. The JC Hamiltonian^i 



+ g{a^(7 +cr+a), 



(Al) 



describes the situation of a two-level system (transition en- 
ergy e, Pauli raising and lowering operators a^) coupled to 
a harmonic oscillator (frequency uj, creation and annihilation 
operators a^, a) with coupling strength g. Prominent realiza- 
tions of Jaynes-Cummings physics include cavity and circuit 
QED^^ where the JC Hamiltonian is usually obtained within 
the rotating wave approximation (RWA) valid for g <C uj,e. 

The RWA is crucial: within this approximation, the total 
number of excitations N = a^a + ct+ct^ is conserved and 
the Hamiltonian separates into 2 x 2-blocks for each fixed- 
excitation subspace, i.e. H"^ = (S)J^o '^m' where ho = and 



Rn 



cos^ 01 sin 01 

i -I i 



^1 + 



(Bl) 
(B2) 



Vn+T cos 9n cos 0n+l + 



In sm Or, sm ( 



Vi+l 



+ 



cos sin 0n+\ \/^ sin O^i cos 0n+l 



I -v/n cos 0n cos 1 + \/n — \ sin 0„ sin 1 1 

(n— 1) — 

|y^cos6'„ sin(?„_i — \/n— lsin0„ cos0„-i| 

^n- ^(n-l) + 



(See the previous AppendixlAlfor the definition of the mixing 
angle 0^.) 



riLj \fng 

\fng (n — l)u; + e 



(A2) 



Precisely because of this block-diagonal form, the JC model 
is exactly solvable, requiring merely the diagonalization of a 
2 X 2-matrix. The resulting energy spectrum is given by the 
ground state energy £'o = and the excited state energies 
(n > 1) 



nu! 



± 



-,1/2 



ng 



(A3) 



Here, we have used A = e — as an abbreviation for the 
detuning, and the polariton number is counted by n G N. The 



APPENDIX C: INVARIANCE OF ATOM-PHOTON 
COUPLING AND SPIN BERRY PHASE UNDER LOCAL U(l) 
GAUGE TRANSFORMATION 

In Section IVll we have claimed that the action 5" given 
in Eq. ( fSSl ) is invariant with respect to the local U(l) gauge 
transformation specified in Eq. ( l63T l. The invariance is verified 
by simple inspection. For completeness, we briefly discuss the 
proof for the atom-photon coupling and the spin Berry phase 
term. 

First, consider the coupling term 



9io.*[N,,^,~iN,,y]+c.c.) 



(CI) 
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This is transformed into 



g [a^e-'"^ [N^, - z7Vj_ J + c.c.^ = g (^a^e"'^ [cos (piV,,, 
+ sin ^pNj y + i sin (pNj^^ — i cos ipNj^y] + c.c.^ 



which reveals the invariance. Second, consider the transfor- 
mation of 5b. Writing the rotation matrix affecting Nj as 
R(r), we have 

Sb^Sb + J2 rdr(N,(r)|Rt(r)R(r)|N,(r)). (C3) 

„■ "'0 



Evaluating the matrix derivative and multiplying with the 
transpose of R, we find 



Rt(r)R(r) = 




(C4) 



which gives (Nj(T)|R^(r)R(r)|N_, (T)) = 0, and thus proves 
the invariance of Sb- 



APPENDIX D: DERIVATION OF RELATIONS BETWEEN 

A'a, ANDf 

Our derivation of the interrelations between the coefficients 
Ki, K2, and f which enter the effective action for the auxil- 
iary field 1}], Eqs. ( |56] | and ( fSTl i. is analogous to the reasoning 
applied in the case of the Bose-Hubbard model; see, e.g., Refs. 
HifandU^ Starting point of this consideration is the invari- 
ance of the action S", Eq. ( ISSl l, under the generalized gauge 
transformation specified in Eq. ( l63T l. Once the Bose fields a^, 
a*, and the spins Nj have been integrated out, this invariance 
remains intact in the following sense: when transforming the 
auxiliary fields and the cavity frequency according to 



{uj-^)^{uj-fi)-i^, (Dl) 
ar 



the effective action S' with Lagrangian density 



remains unchanged. This statement must hold for arbitrary 
configurations of the fields ip* and tpj, as well as arbitrary 
choices of the phase function ip{t). As a consequence, the 
invariance must not only hold for the effective action, but also 
for the Lagrangian density itself. In particular, the coefficients 
Ki, K2 etc. themselves must be invariant, and they must be 
so to each order in Lp, ip, and its higher time derivatives, when 
considering a small global phase change (p{t) <C 1. 

Carrying out this expansion for the terms included in £eff 
above, one finds 



£eff - i^K'a - -^^K^ 



(D3) 



2 



lip- 



dtp 

iipK2 ) V'*-^ + c.c. 



(A'2 - i^k:^) 



dip 



dr 



+ (r — iipr + iKiip) + —{u — i(pu') j?/'! + 



where Cp ^ d^j dr, and primes denote derivatives with respect 
to {u!—ji). It is important to note that the terms collected in Eq. 
( ID3l l do not yet correspond to a consistent expansion. Specif- 
ically, higher-order terms dropped in £eff above give relevant 
contributions even to the order recorded in Eq. ( ID3l l. As one 
example, consider the higher-order term 



dip* d'^ip 
dr dr'^ ' 



(D4) 



whose expansion contains (among others!) the contribution 



dip 



(D5) 



One can verify, however, that the coefficients of the \ip\ and 
the ip*dip/dT terms are in fact consistent to the order 0{ip) 
as written above, and we can read off the important relations 



1 f ^dip dip* ' 



) +K2 


dip 




dr 



(D2) 



(D6) 
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